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Quantum error correction (QEC) |1] is fundamental for quantum information processing but entails a substantial over- 
head of classically-controlled quantum operations, which can be architecturally cumbersome to accommodate. Here 
we discuss a novel approach to designing elementary QEC memory cells, in which all control operations are performed 
autonomously by an embedded optical feedback loop. Our approach is natural for nanophotonic implementations in 
which each qubit can be coupled to its own optical resonator, and our design for a memory cell based on the quantum 
bit-flip or phase-flip code requires only five qubit-cavities (three for the register and two for the controller) connected 
by wave-guides. The photonic QEC circuit is entirely on-chip, requiring no external clocking or control, and during 
steady-state operation would only need to be powered by the injection of constant-amplitude coherent fields. 

While our approach exploits rigorous theoretical results derived using techniques from quantum field theory and 
coherent-feedback quantum control [2H7], we devote the main body of this article to an essentially intuitive presentation 
of our QEC photonic network. Following this, a more complete modeling of the system is presented in extended 
appendices. Appendix [A] briefly recounts the Hudson-Parthasarathy stochastic calculus [3H5] and a modeling of 
interconnected open Markov quantum components [5] that provide a formalism well-suited to modeling complex 
quantum optical networks. Appendix|B]then provides the specific models for several of the components in our network 
and introduces a key limit theorem [3] necessary to produce a computationally-tractable model for the overall network. 
Finally, Appendix [C] constructs the full network model and derives its dynamics. 

The photonic circuit shown in Fig. [l] implements continuous QEC based on the bit-flip or phase-flip code 1 . Ql, 
Q2 and Q3 represent the register qubits; the blue signal lines indicate the routing of a laser beam (injected in the 
coherent state |2a), where |a| 2 has units photons time -1 ) used for error detection (following principles described 
in [5]); i?l and R2 represent two qubit-based photonic relays (as analyzed in [5]); and the red signal lines indicate 




FIG. 1: Photonic QEC network schematic. Schematic diagram of a coherent-feedback quantum memory showing qubits 
in cavities (Ql, Q2 and Q3), circulators, beam-splitters, steering mirrors and relays (Rl and R2). 
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FIG. 2: The cQED qubit. Level diagrams of a, Z probe interaction, b, X probe interaction. 
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FIG. 3: The cQED photonic relay. Details of the cQED relay component model (modified from [9]): a, input and output 
ports, b, coupling of input/output fields to resonant modes of two cavities, and c, relay internal level diagram. 



pathways for optical feedback using two additional laser beams (injected in the coherent state |/3)) for corrective 
Raman bit- or phase-flips. Such a circuit could be implemented using solid-state qubits embedded in nanophotonic 
resonators connected by single-mode wave-guides |101 lllj , or using microwave components as in circuit QED |12) . 

The quantitative performance of such a quantum memory will of course depend on physical parameters such as the 
qubit-cavity coupling strengths and the cavity decay rates, but below and in Appendix [B] we show that our approach 
does not require any fine-tuning of the parameters and requires only that certain ratios be large. It is interesting 
to note that these ratios become large in a small volume limit for the optical resonators, which is a natural limit to 
consider in nanophotonic implementations. Passage to the small volume limit also gives rise to drastic simplifications 
of the quantum input-output models for the components in the QEC circuit, and thus emerges as a significant new 
abstraction principle for the analysis of nanophotonic signal-processing networks. 

The physical qubits Ql, Q2 and Q3 are actually multi-level 'atoms' or comparable solid-state emitters, e.g., nitrogen- 
vacancy centers. Fig. [2] shows the level structure presumed for the individual atoms and the couplings utilized for the 
bit /phase-flip schemes, with {\g), \h)} corresponding to the single qubit states {|0), |1)}. In the bit-flip implementation 
we assume that each qubit is strongly coupled to a single-mode, single-sided cavity on its \g) <-» \e) transition as 
discussed in [5]. As a result, the phase acquired on reflection by a probe beam is 7r if the qubit is in state \g) and if 
it is in \h). While details of the probe-qubit interaction can be investigated for any desired physical parameters using 
the familiar cavity QED master equation [8], for our network calculation we utilize a simple abstraction obtained 
when the vacuum Rabi frequency g p and cavity field decay rate k p are taken to infinity with the decoherence rates 
held fixed. Because g p and k p both increase with decreasing resonator volume, we term this a 'small volume' limit. 
The simplified model is derived using a quantu m sto chastic differential equation (QSDE) limit theorem [13. following 
a similar procedure to that of [5J (see Appendix B 1 1 . In the phase-flip scenario we modify the scheme by assuming a 
linearly polarized cavity mode (as in a micro-disk or photonic crystal resonator) that couples to both a cr_ transition 
from \g) f-> |e) and a a+ transition from \h) «-» |e). 

The central component in the network's QEC 'controller' is a recently proposed [9] cavity-QED relay, some details 
of which are presented in Fig. [3j Three cavity modes are resonant with various transitions in an intra-cavity atom, 
as shown in Fig. |3]b and c. In a similar small volume limit [5], the atomic dynamics are limited to only the ground 
states, and when the atom is in state \h) (\g)) the cavity's POWER input is connected to the OUT (OUT) output 
port. A probe signal driving the SET input port (in the absence of signal at the RESET input) causes the relay state 
to decay to \h), while driving only the RESET input induces decay to \g). Through these internal dynamics, a probe 
signal entering either the SET or RESET port of the relay routes the POWER input out either the OUT or OUT 
port, respectively. 

We now trace through the circuit dynamics of Fig. [T] in more detail, assuming a bit-flip code implementation. A 
coherent input \2a) is split into probe and reference beams. The probe beam first reflects from the cavity containing 
Q2, gaining a ir or phase shift depending on the state of that atom, putting the probe signal into a coherent 
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state with amplitude \/2aZ^ 2 ^ (where Z^ 2 ^ = \h){h\ — \g)(g\ on Q2 and the identity operation on all other Hilbert 
spaces is implied — see Appendix [C]). The probe signal is again split and the two resulting fields separately probe the 
cavities containing Ql and Q3, resulting in signals with coherent amplitudes qZW 1 '^ 2 ' and aZ^ 2 ^ Z^ 3 \ These 
are interfered with copies of the reference beam, producing four signals with amplitudes a(I ± Z^^Z^ 2 ^)/^ and 
a(I± Z^Z^)/y/2. The '+' signals are sent to the SET inputs of relays Rl and R2, the '-' to the RESET inputs. 
In this way, the SET (RESET) input to Rl receives a coherent input signal if Ql and Q2 have even (odd) parity, while 
the RESET (SET) receives only a vacuum input. Similarly for R2 and qubits Q2 and Q3. The POWER inputs to Rl 
and R2 are two more coherent fields, each with amplitude f3 and frequency and polarization that drive the remaining, 
far-detuned atomic transitions depicted in Fig. [5J The OUT port of Rl (R2) is directed to both Ql and Q2 (Q2 and 
Q3). The OUT port of Rl (R2) is directed to Q3 (Ql). When a qubit is simultaneously illuminated by two of these 
feedback beams, Raman resonance conditions are satisfied, inducing a coherent Rabi oscillation between the states of 
that single qubit. When only one feedback beam illuminates an atom, the beam is assumed to be sufficiently detuned 
such that it only imparts an AC Stark shift to one of the atomic ground states. Remarkably, a proper choice of the 
ground state each feedback beam targets prevents Stark shift-induced dephasing of the stored codeword in any error 
state (see Appendix C 4 ) . 

Assuming the register qubits are initialized in a proper bit-/phase-flip codeword, this network will self-stabilize 
against bit /phase-flip errors without any external input other than the coherent probe and Raman fields (a, j3, and 
all Hamiltonian coupling constants are constant in time). For an intuitive understanding of the network, note that 
if no flip errors occur, the SET inputs to the relays are activated, Ql and Q3 only receive one feedback signal each 
and no corrective dynamics occur. If for example Q2 flips, the RESET inputs on both Rl and R2 are activated, Q2 
receives both feedback signals, while Ql and Q3 again get only one each. A corrective flipping Hamiltonian is thus 
implemented on Q2 until a flip is 'complete,' which in turn drives the relay switches back to the 'no-error' state, which 
then extinguishes the flipping Hamiltonian. This coherent feedback network thus implements a continuous version 
of the usual discrete-step strategy that diagnoses and corrects particular errors without disturbing the information 
content of the encoded joint state. 

After interconnecting simplified (small-volume) models for the bit-flip components according to the bit-flip circuit 
diagram, including individual bit-flip error processes, and adiabatically eliminating of all but the qubit and relay 
ground states, we obtain (see Appendix [Cl the closed-loop master equation 



p t = -i[H,p t ] 
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where X^ = \h)(g\ + \g)(h\, — and 11^- = |j)(j| on component C, each qubit suffers random flips 

at a rate T, and the 'feedback parameter' is proportional to |/3| 2 /A (note that several Lj terms corresponding 
to strong dephasing of the relay states have been omitted as they have no appreciable effect on our dynamics). 
An analogous procedure using the phase-flip component models leads to a master equation with the substitutions 
ZiQi) o XW), 2ni Ql) ->■ 1 ± iW'l While L 5 through L 7 drive the corrupting bit-flip processes, the first three 
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terms of the Hamiltonian (H) and remaining 'collapse operators' work in concert to implement the corrections: L\ 
through L 4 chiefly drive the state of the two relays to reflect the qubit parities, and each bit-flipping term in the 
Hamiltonian is only active for a particular 'syndrome state' [TJ, as indicated by the relays. The remaining Hamiltonian 
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FIG. 4: Network performance 

(see Appendix C3bl, a = Q/8 
is also shown. 



Decay of fidelity, {^o\pt l^o), for several values of the feedback parameter Q = |/3| 2 7/2A 
and r = 0.1. The fidelity decay of a single, bare qubit also suffering bitflip errors at rate r = .1 



terms physically correspond to AC Stark shifts, whose effects have been minimized by design (see Appendix C 4 1 . 
Indeed from the general structure of equation |2]) it is straightforward to intuit master equations that correspond 
abstractly to memories based on other stabilizer codes, highlighting the potential of our approach for developing 
technology-inspired abstractions in engineering of quantum optical networks. 

Note that a network simulation with full physical models of all the components would be beyond our present com- 
putational capabilities; judicious use of the small- volume abstraction, which allowed us to derive the easily interpreted 
closed-loop master equation ([T]), has thus been essential in our effort to prove that autonomous QEC with fully em- 
bedded control is possible in principle in nanophotonic implementations. Previous studies [HI [5] have shown that the 
idealized behaviors are reasonably well approximated in component simulations using realistic physical parameters. 

In Fig. [4] we display some illustrative numerical integrations [14] of equation ([!]) for the bit-flip scenario. With 
the flip rate set to Y = 0.1, the feedback amplitude is varied from O = (no feedback) to 57 = 210, with the 
probe amplitude kept at a — The initial state of the qubit register is chosen as |\&o) = (\ggg) — i\hhh))/y/2, 

Pa = | v fo)(*oli an d the fidelity decay (^olPtl^o) is computed to quantify the decoherence suppression achieved by 
the coherent feedback loop. With f2 — the fidelity decay is identical to that of the three-qubit register without any 
measurement or feedback, while for larger values of fl and longer storage times, the fidelity of the encoded qubit is 
clearly improved relative to what it would be for a bare qubit. Variation in performance for different l^o) is only 
due to the different susceptibilities of codewords to bit-flip errors. Because of the small, but finite likelihood that 
additional errors accrue before the correction of a single flip is implemented, fidelity is steadily lost for finite probe 
and feedback strengths, in analogy to the failure of the discrete time bit-flip code to correct multiple flips. 

In conclusion, we have described a novel approach to continuous quantum error correction based on stationary 
coherent feedback. Our scheme may match well with hardware implementations that feature solid-state qubits em- 
bedded in planar circuits of electromagnetic resonators and wave-guides, and establishes important new connections 
between QEC and coherent-feedback quantum control theory [1514T7] . While our current results are limited to the 
bit /phase-flip scenarios, we hope in future work to find ways to combine bit-flip and phase-flip probe interactions 
(possibly employing either the five-qubit code [T] or Bacon-Shor subsystem code |T5] ) to design feedback networks 
that can correct arbitrary single-qubit errors. 
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Appendix A: QSDE formalism 



To describe our QEC coherent-feedback network, we employ the formalism of the Hudson-Parthasarathy stochastic 
calculus and quantum stochastic differential equations (QSDEs), and the concatenation and series product rules for 
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interconnecting open Markov quantum components. Here we shall only give a brief reminder of these formalisms, a 
detailed exposition and nice introduction to the formalisms, including several illustrative examples, can be found in 
[HE]; see also the references cited therein. 



1. Hudson-Parthasarathy QSDEs and unitary evolution of open Markov quantum systems 



In the following we shall use the notations: 



1, * will denote the adjoint of a linear operator as well as the 



conjugate of a complex number, if A — [ajk] is a matrix of linear operators or complex numbers then 
and A> is defined as A> = (A#) T , where T denotes matrix transposition. We also define Re{A} = (A + A#)/2 and 
Im{A} = (A — A#)/2i and denote the identity matrix by / whenever its size can be inferred from context and use /„ 
to denote an n x n identity matrix. 

Let r]l,...,i]™ be n independent vacuum quantum white noise processes satisfying the commutation relations: 
Wt > (.Vt 1 )*} = $jk5(t — f) and \rf t ,r) k ,] — Vj, k and Vt, t' > 0, and define A{ — L rf s ds (j = 1, . . . , n) to be a vacuum 
annihilation field, A{* = J *(?7s)*^ s to be a vacuum creation field, and A{ k = J* (r] 3 s )* r] k ds . The processes A\, A{* , A\ k 
are quantum stochastic processes collectively referred to as fundamental processes, and A\* = (A\)* . The products of 
their differential increments are given by the quantum Ito table 01 [5] 
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All products dMt dt and dt dM t vanish. 
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If two processes X t and Y t have the quantum stochastic differentials: 



dX t = 
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dY t = 



F 2j dAj + G 2j dA{ + J2 H 2jk dAf 



(Al) 



k=l 



for constant system operators Fij,Gij, Hijk, I = 1,2, then the quantum stochastic differential of the product X t Y t is 
given by the quantum ltd rule: 



dX t Y t = dX t Y t + X t dY t + dX t dY t , 



(A2) 



in which each term on the right hand side can be evaluated using the quantum Ito table given above, and using the 
fact that X t and Y t commute with each of the differentials dA 3 t , dA{* and dA{ k [H|S]. 

Let A t = (Al, A\, . . . , A") T . In the interaction picture with respect to the free field dynamics, the general evolution 
of an open Markov quantum system coupled to the fundamental processes at any time t > is given by a unitary 
(co-cycle) propagator V(t) satisfying the right Hudson-Parthasarathy (HP) quantum stochastic differential equation 
(QSDE) @HB]: 



dVt = 




rjk 



S jk )dA{ k + dA\h r - L\SdA t - (iH r + -L\L r )dt V t ; V(0) = I, 



(A3) 



is self-adjoint (i.e., Hermitian) linear operator that is the Hamiltonian of the system, L r = 



where H,, 

(L ri i,i r 2, ■ • • ,L rt7l ) T is a column vector of coupling linear operators L r j to the fields, and SVjfc, j, k = 1, . . . ,m, 
is a bounded system operator such that S r = [S r ,jk]j,k=i,...,n is unitary, that is, S$S r = S r Sl = I. S r is called a 
scattering matrix. Here we have added the subscript r just to emphasize that (A3) is a right HP QSDE. If we define 



K r = -iH r - \ Y^k=\ L r k L r,k, then we see that K T + K* = - J2k=i L *r k L r,k and K r - K* = -2iH r . Using the right 
propagator V(t), observables evolve in the Hciscnbcrg picture according to a t = V t * a$Vt, which has the differential 
form 



da t = V* a dV t + dV t * a V t + dV t * a dV t . 
Note that the right QSDE is completely specified by three parameters, S r , L r , and H r 



(A4) 
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In many parts of the calculations to follow in these appendices we will work with left QSDE's involving a left 
unitary propagator U t such that observables evolve according to a t = U t aoU t * . We use this left form because we will 
be applying the adiabatic elimination results of Bouten, van Handel and Silberfarb 13J, in which they adopted the 
left QSDEs in their work for technical reasons, rather than physical ones; see [T31 Remark 3]. 

The relation between the left propagator Ut and right propagator Vt for the same system is 

dV t =dU?. (A5) 



2. The master equation 



Working with the right QSDE we follow Gardiner's method (see, e.g., [7j) for deriving the master equation; see also 
the derivation in [U |S] . By averaging over the noise terms in the evolution, we obtain the expected evolution of an 
observable 



d (a t ) = (dot) = (V t * a dV t + dV t *a Q V t + dV* a dV t ) 
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(A6) 
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(A8) 



3. Concatenation and series product 



Let G denote an open Markov quantum system that evolves according to the right QSDE (A3 ) with given parameters 
S r , L r and H r . For compactness, we shall use a shorthand notation of [2 , and denote such a system as G = 
(S r , L r , H r ). In this section we briefly recall the concatenation and series product developed in j^j that allows one to 
systematically obtain the parameters of an open Markov quantum system built up from certain interconnections of 
open Markov quantum systems. 

Let G\ = (Si, Li, Hi) and G2 — (5*2,^2)^2) be two open Markov quantum systems. The concatenation product 
G\ EH G2 of Gi and G2 is defined as 



Gi EH G-2 — ( 



'Si 0" 




' Li ' 


s 2 _ 




. L 2. 



,Hi+H 2 ). 



(A9) 



It is important to note here that elements of Si,Li,H\ need not commute with those of <S2, £2,-^2 ^' e "' ^ 1 an< ^ 
may be the sub-systems of a common larger system). If Gi and G 2 are independent systems, in the sense that their 
coefficients commute with one another then concatenation can be interpreted simply as the grouping of variables of 
two non-interacting open Markov systems to form a larger open Markov system. 

It is also possible to feed the output of a system Gi (here the output is the field that results after the incoming 
field has interacted with the system; see [H [7] for details) to the input of system G 2 , with the proviso that Gi and 
G2 have the same number of input and output channels. This operation of cascading or loading of G2 onto Gi is 
represented by the series product G2 < G± defined by: 

G 2 <iG 1 = (S 2 Si,S 2 L 1 +L 2l H 1 +H 2 + lm{LlS 2 L 1 }). (A10) 

Note that G2<Gi is again an open Markov quantum system with a scattering matrix, coupling operator and Hamilto- 
nian as given by the above formula. Moreover, there are two important decompositions of an open Markov quantum 
system based on the series product. They are 

(S,L,H) = (I, L, H) < (S, 0, 0) = (5,0, 0)<(/,5 t L, H), (All) 
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4. Coherent fields and the Weyl operator 



Thus far, we have provided overviews for QSDEs driven by quantum noise fields that are in the vacuum state. 
However, in many situations in quantum optics the fields are in coherent states, emanating from a laser source. 
Suppose that each of the fields A], A 2 , ... , A™ is in a coherent state with (complex) amplitude at, cx.2, ■ ■ ■ , otn, and 
let a = {at, 012, ■ ■ ■ , ct n ) T ; note that this includes the vacuum states by setting the amplitude for a particular field to 
zero. Fields with coherent states can be obtained from fields in the vacuum state by an application of a Weyl operator 
satisfying both the right and left QSDE: 

dW t a = (dA\a - aUA t - ha\ 2 dt)W t a , 

= W t a (dA\a - aUA t - \\a\ 2 dt), (A12) 

and we also have the relation 



W t a * = Wf a . (A13) 

Then the right unitary propagator V t a of a system driven by fields with coherent amplitudes a is V t a — VtW", 
where Vt is the right QSDE given in (A3), representing the same system with vacuum inputs. Alternatively, using 
the series product formalism from the preceding subsection, we can write G a = G < (I, a,0), where G a — {S,L r + 
Sa, H r + Im{L' Sa}) is the system driven by coherent fields with amplitudes a, and G is the system when driven by 
vacuum fields. 



Appendix B: Component and subsystem models 

In this section we describe several of the network subsystems using the formalism introduced in section [XJ We also 
introduce and apply a recent adiabatic elimination theorem |13j that generates a QSDE description of an often much 
simpler, approximate evolution of a physical system with well-separated time scales. These approximate models form 
the subsystem components of the full network model presented in section [Cj 



1. Probe component model 



In this section we present the physical cQED model for the probed qubits and derive its limiting evolution used in 
the network model [8]. 

As represented by the blue transitions in Fig. [2j the Z and X probe interactions emerge from a three level atomic 
system in a A-configuration interacting with a quantized mode of a single-sided electromagnetic resonator (cavity) . In 
the case of the Z probe, the \g) O |e) transition is strongly coupled and on resonance with the cavity mode, while the 
\h) state is fully uncoupled. In the case of the X probe, both \g) o \e) and \h) |e) transitions are simultaneously 
coupled and on resonance with the cavity mode. The X probe configuration could be realized with a degenerate 
spin-i ground state pair coupled to an excited spin-0 state via a linearly polarized cavity mode, for example. Thus, 
in a properly rotating frame and applying the usual rotating wave approximation, the familiar Jaynes-Cummings 
Hamiltonian interaction for both systems is 

H = ig c (a*a — era*) (Bl) 

where a is the cavity mode annihilation operator, the atomic 'lowering operators' act as a — |<?)(e| on the atomic 
states in the bit-flip network, a = (\g) + \h)) (e\ in the phase-flip network, and the coupling rate g c may be taken 
real without loss of generality. 

However, the cavity mode and atom also interact with the free fields that surround them: the cavity photons will 
decay out of the resonator and the excited atomic state will spontaneously emit photons into the environment. These 
interactions are included in a HP QSDE description of the dynamics with the inclusion of two coupling operators [19] : 

Li = \/2na] L 2 = y^r/Iv (B2) 

where k and j± are the cavity field and atomic dipole decay rates, respectively. As there is no direct coupling 
between the cavity and atomic decay channels, the scattering matrix for this system is the identity. Thus, the right 
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QSDE descriptions of the dynamics for the Z and X systems with vacuum inputs are characterized by the operator 
coefficients 



Q 



2na 



,tg c [a a- 



aa 



(B3) 



We will eventually use the limiting evolution of these 'vacuum' cQED systems in our network model, but note that 
familiar, driven cQED dynamics may be derived with the inclusion of a Weyl operator in series product with the 
vacuum system (representing an on-resonance, CW coherent-state input to the cavity mode) 



,0 



2na + a 
27-1 cr 



, ig c (a*a — era*) H — (aa — aa*) 



(B4) 



followed by the application of the master equation from section A 2 



We now begin an adiabatic elimination procedure on Q to obtain the very simple limiting dynamics of a measurable 
qubit, which will be used in the network model. The physical intuition of the derived dynamics is as follows. As 
described in [5], in a bad cavity limit (n — > oo), with the atom in a uncoupled state, a coherent probe resonant with 
the cavity will enter the cavity and exit 'immediately' (due to the large k), picking up no additional phase shift. 
However, if the atom is in a coupled state, and the system is also in a 'small volume' coupling limit (g c , k — > oo, 
with g c / k constant), then the same probe will simply reflect from the off-resonant, atom-cavity system, picking up an 
additional ir phase shift. Thus, in these limits, the cQED system acts as a simple scattering object with no internal 
dynamics, but which imparts a phase shift on a cavity-incident coherent probe, as determined by an internal qubit 
state. 

Following the mathematical formalism of [13] . we redefine and relabel Q with a scaling parameter k, which is 
assumed to take the limit k — > oo, and write the component as a left QSDE 



dQ ( t k) = Q^j - kV2^adAl* + kV2na*dAl - -j2^adA 2 * + v / 2^a*dA 2 - 

k 2 na* adt — j±a*adt — k 2 g c (a* a — aa*)dt^. 

Note that with the scaling, the cavity field decay and coupling rates take the limits described above. 
In the notation of |13j . the operator coefficients of the above left QSDE are 



= -k 2 na*a - - k 2 g c {a* a - acr*) , L {k) 



k\l2na* 



(B5) 



Through the adiabatic elimination procedure of [13] , a set of limiting operator coefficients (which we denote without 

the '(*>') are derived, which define another, limiting propagator, Q t , to which strongly converges. It is this 
limiting propagator that will be used in the network model in place of the physical cQED model. 
For the adiabatic elimination [T3J section 2.2] we use 



Y = —na*a — g c (a*a — aa* 



F : 



2na* 







G = 







/27la* 



A = 0, B = -7_LO-V, 

W = I. 



(B6) 



In the Z probe case for the bit-flip network, we choose Hq = span {| g, 0) , | h, 0)} and define 



Y{\g0),\h0)} 


= 0, 




Y \ hn) 


= -—\hn) 
nn 


j n>l; 


Y \e (n-l)> 


_ 9c- 


Jn 


K 2 n (n — 


1) + gin 


Y \gn) 


k (n 


-1) 


K 2 n (n — 


1) + gin 



g 2 c n 



nn 

K 2 n (n — 1) - 
g c \fn 



K 2 n (n — 1) + g 2 n 



\e(n-l)), n>l; 



|e(n-l)), n>l. 



(B7) 
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This provides YY = YY = P±. Then by [13], the limiting operator coefficients are, 

K = 0, Li = L 2 = 0, Mi = M 2 = 0, 

and (with Tl in = |i,n)(«,n|) 

AT n = p Wu (f*YF 1 + <5 U ) P = P + 2 K P af a* (n g0 + II M ) 
= P + 2 K P oa y (\gl) (g0\ + \hl) (h0\) 
= P + 2«P a [ — I eO) (g0\ --\hl) (h0\ 

\9c K 

= P -2\hG) (h0\ 

= n s o — ii/jo = z, 



(B8) 



N 12 = P Wu (F*YF 2 +S 12 )P =0 



n 21 = p w 22 f:yf 1 + s 21 p = o 



n 22 = p w 22 [f;yf 2 + 5 22 j P = 1, 

Collecting the results, the limiting system is characterized by the operator coefficients 



Q = 



z o 
o I 



,0 



(B9) 



(BIO) 



to which the physical cQED system, Q^ k \ strongly converges in the limit k — > oo. For the X probe case in the 
phase-flip network, the analogous follows, but with the pauli-X spin operator eventually replacing Z in equation 

pTob. 



Raman transition model 



In this section we describe the physical Raman interaction that is used to correct the qubit states in the coherent 
feedback loop. We presume in this case that the atoms in the probed qubits described in subsection Bl are also 
interrogated by two free field modes (i.e. field modes transverse to the cavity mode) that drive atomic Raman 
transitions when displaced by a large, coherent amplitude. 

We assume that there also exists an additional atomic excited state |r), distinct and far detuned from \e) for 
theoretical simplicity, although \r) = \e) should work in practice. Both ground states are coupled to |r) via interactions 
with transverse free fields. Two of these field modes, far detuned by frequency A from the {\h), \g)} O \r) transition 
are used to drive the Raman interaction, as represented by the red transitions in Fig. [2j For the bit-flip network, 
in a properly rotating frame, these interactions are modeled in a QSDE by the Hamiltonian and coupling operators 
(redefining A — !• k 2 A in anticipation of an eventual k — > oo limit) [20] 

H {k) = fc 2 A|r)(r| = /c 2 An r 

4 k) = yfl\h){r\=^c7 hr 

L 2 k) = Vl\9)(r\ = ^a gr . 

For the phase-flip network, the above coupling operators are modified with Ohr — > ~^ 2 \ (J hr + fgr) arid <j g 
<j gr ) (physically suggesting two linearly polarized modes coupling to angular momentum eigenstates, for example). 
Also, the interactions of the {\h}, \g)} -f-> |r) transitions with all other free modes (which principally drive spontaneous 
emission decay from \r) to the ground states) may be aggregated into just two more coupling operators 



(Bll) 



75 ( a hr 



V2 



rW _ 
^3 - 

r(fe) _ 



y\\<7hr 



(B12) 
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in the limit 7 <C 711 ■ Together with an identity scattering matrix, the above Hamiltonian and coupling operators 

represent a simple vacuum- input Raman system, . 

In the quantum network, the field modes 1 and 2 undergo a non-trivial processing before reaching the Raman 
systems described above. For demonstration purposes, however, consider a bit-flip Raman system in isolation, but 
with modes 1 and 2 simply displaced into large-amplitude coherent states. This system can be represented by the 
operator coefficients 



y8,(*) 



<(w m mw m m{i 2 ,o,o)) 



( 



1, 



V 



y/l<Tgr + k(i 2 



, fc 2 An r 



iky/7 



\ 



{ogrfil - cr* r /3 2 ) 



(B13) 



/ 



where W k p is a Weyl operator of displacement amplitude kfi. 

Application of the adiabatic elimination theorem |13] to the above (using a technique discussed below in section 



B 3 ) results in an approximate Hamiltonian interaction 



H = --jr{\h\ 2 U H 



l/3 2 | 2 n 9 



(B14) 



and all atomic operators drop out of the field coupling operators (eliminating all spontaneous emission dynamics). 
Thus, coherent transitions between the two ground states may be driven by the last two terms in H when both Raman 
modes are in large-amplitude coherent states. However, the ground state energy shifts represented by the first two 
terms (physically recognizable as AC Stark shifts) are generally inconvenient in the QEC network. 

To compensate for these Stark shifts, in principle we could presume an additional pair of excited states {\H}, \ G)} 
(not shown in Fig. [2| that also couple to the ground states {\h}, \g)} via modes 1 and 2, respectively, but with the 
opposite detunings: — fcA. This choice to add more atomic levels is again made for theoretical simplicity, although 
two additional free field modes coupling to |r), but with opposite detunings, would likely be more practical. Thus, 
the physical Raman interaction would have the form 



Qr+ — ( ^2; 



27(07^ 
27(cr S r 



<?Hr) 



,fc 2 A(n r -n H -n c 



(B15) 



where we have omitted the coupling of the excited states to any other free modes, as the dynamics they induce would 
eventually drop out in the k — > 00 limit, as above. In this limit, the effective evolution is again governed by the 



Hamiltonian dynamics depicted in equation (B14|, but without the first two, AC Stark shift terms. 



3. Adiabatic elimination with large held displacements 



In the previous section we mentioned that the adiabatic elimination processes will be applied to the Raman com- 
ponents only after they are placed in the quantum network. However, there are some mathematical complications in 
applying the usual theorem [T3] when the amplitude of coherent field input must be scaled. In this section we describe 
a procedure to overcome this difficulty briefly and in the abstract. This technique is a generalization upon a similar 
technique developed by Luc Bouten and applied in [8]. 

Let (S° , L° , H ) be a network with only vacuum inputs. The same network, but with coherent displaced inputs is 
defined as 

G« = (S Q ,L ,ff )« (7,^,0) 

= (S°,S d k +L ,H° + Q{L^S°d k }) = (S,L^,H^) (B16) 

where is a vector with complex- valued entries, some of which may scale with k (e.g. kfl), some may not (e.g. 
2a), and some still may be 0. The scattering matrix S = S° is assumed to have no scaling with k. It happens that 
applying the adiabatic elimination theorem on is complicated by the S°dk term in the coupling operators. The 
first step in overcoming these difficulties is to reconstruct the network such that all the scaled parameters are moved 
to the Hamiltonian of a new subcomponent. That is, we rewrite 



G (fe) = (S, 0,0) < (I, d k ,0)< {I, S^L°,H^ -3{4s ot L }) 
= (5 ! 0,0)<(/,4,0)<GW 



(B17) 
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It will be straightforward to apply the adiabatic elimination procedure on the subsystem & k >, which will produce a 
limiting subsystem G. Thus G^ is asympotically approximated by 

G (fe) = (S, 0,0) < (l,d k ,0) <G (B18) 

for sufficiently large values of k. It is this asymptotically approximate and simpler network & k ' that will be used in 
final description of the coherent feedback system. 



Appendix C: The quantum network 



Here we present a detailed calculation of the coherent-feedback network model. A schematic diagram of the overall 
quantum network is shown in Fig. [5} We start by introducing notations and reviewing the component models; detailed 
derivations for several of the component models are collected in section FBI 
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FIG. 5: The quantum 'correction' network model. Concatenation decomposition of a quantum network analogous to 
that represented in Fig. [lias G EB G' . See Remark [l] for an explanation of the network differences. 



Remark 1 The network considered in the appendices is slightly different from the network shown in Fig. [7j This 
is because in these appendices, the cavity for Q2 is taken to be two sided with a coherent probe laser impinging on 
each side, allowing us to obtain a decomposition of the network into two (symmetric) parts as shown in Fig. [5] and 
simplifying the calculations for the overall network parameters that are presented in subsection C3 On the other 
hand, in Fig. [7] we make all qubit cavities identical and one sided, which would be probably more suited in actual 
cQED implementations. Similar effects can be achieved with either a double sided cavity for Q2 or a single sided 
cavity that is followed by a 50/50 beamsplitter. Although the two networks are slightly different, they only differ in 
the overall network scattering matrix, and a slightly modified repeat of the lengthy calculations of subsection \ C ffi with 
Q2 one-sided and followed by a beamsplitter, will show that the network coupling vector and Hamiltonian will be the 
same for both networks. As a consequence, we emphasize that the two networks will have the same network master 
equation. Since, our main interest here is the network master equation, to derive this equation it is more convenient 
to work with the slightly modified network of Fig. [5[ 
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1. Notation 



We shall use the following symbolic convention for individual components in the network: B = beam splitter, Q = 
qubit cavity, R = quantum relay/switch. Distinct network components that are of the same type are distinguished 
by assigning to each a unique subscript as an index, and all components in the network are open Markov quantum 
systems described by three parameters introduce in section [X| a unitary scattering matrix, a coupling vector, and 
Hamiltonian operator. Matrices and vectors can have operator entries and an entry of a matrix M in row k and 
column I is denoted as Mm- A network component Zy. with Z £ {B,C,R} and k some index will be written as 
Zk = {S^ Zk \ L^ Zk \ i?( Zfc )), where S^ Zk \ L^ Zk \ H^ Zk) are, respectively, the the scattering matrix, coupling vector, and 
Hamiltonian operator of Zk- 



If the matrix S {Zk) is diagonal of the form diag(sjf fc) , S. (Zfc) 



q(Zk)\ 

■ , Onn j 



(here diag(Mi,M 2 , . . .) denotes a block 



diagonal matrix with block diagonal elements Mi, M2, 



with Zi. 



(S( z «), #(^')) and 



S 



(Zk) 



. .) then Zk can be decomposed as Zk = Zj-i EH Zki EH ... EH Zk n 

, L^i) = hf k) and H^ Zk ^ satisfying ELi# (Zfej) = H^ Zk \ 



Jkj v ~ , - , — j -~- ~ -^ , - ~ ~ „„™^.„ 

We will exploit such a concatenation decomposition whenever they can be made to simplify the network calculations. 
Recall that it is not required that entries of the parameters of Zkj commute whenever j ^ k and series products can 
be formed from such components. 

It is important to note that the series product is only defined for systems that have an equal number of input 
and output fields. However, naturally it is often necessary to make interconnections between systems with different 
numbers of input and output field fields. To apply the series product we augment the systems to be connected with 
"dummy" trivial couplings to additional fields as required to make the number of input-output fields of systems to 
be connected be the same. This is best illustrated with a simple example. Suppose we wish to connect a system 
Ci = (S^ Cl \ L (Cl \ H^ 1 ^) having only single input and output fields with another system C 2 = (S (C2) , L^ C2 \ H^) 
having two input-output fields by passing the output field of C\ as the second field into C%. Then two systems 
connected in this way is given by the series product C 2 < ((1, 0, 0) EH Ci)), where (1, 0, 0) represents a dummy coupling 
of G\ to a second field that merely "passes through" C\, without interacting with it, before entering C 2 as the first 
input field. In general, all that is needed is to insert these trivial (1, 0, 0) systems in the appropriate locations before 
making the series connection; see [5] for several examples. 

Recall that /„ denotes the n x n identity matrix and / and to denote identity and a matrix will all zero entries, 
respectively, whose dimensions can be determined from context. In order to reduce the complexity of expressing 
trivially augmented systems we introduce the additional notation EH^ to denote the operation of augmenting a system 
on the left of EH& by inserting a trivial (J„,0,0) system as the fc-th to (fc + n — l)-th component of the system. To 
illustrate, for the system C2 that we had before we would have: 



(C 2 EH 2 (1,0,0))EH 4 (1,0,0) = 








o(C 2 ) - 
D 12 
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c(C 2 ) 
°21 





o(C 2 ) 
°12 




r(C 2 ) 
L L>2 J 


q(C 3 ) 





q(C 2 ) 

D 12 


0" 
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c(C 2 ) 
°12 







r(C2 
L 2 











1 









(CI) 



and so on in this fashion. 



2. Network components 



a. Qubit cavities 



We note that in our quantum network all qubit cavities will have a concatenation decomposition. In what follows, 
Q2 will be taken to be a double ended cavity which is connected to two probe lasers and two bit correcting lasers (cf. 
Remark [T]). It has a concatenation decomposition: Q 2 = Q21 EH Q 22 EH Q23 EH Q 2 4, where Q21 and Q24 are subsystems 
associated with the coupling of Q2 to the probe lasers, while Q22 and Q23 are associated with the coupling of Q2 to 
the bit correcting lasers that drive qubit Raman transitions. Ql and Q3 are single ended cavities decomposable as 
Qj = Qji EH Qj2 EH Qj^ (j = 1, 3) with Qji and Qj$ associated with the coupling of Qj to the two bit correcting lasers, 
while Qj2 is the subsystem associated with the coupling of the cavity to a probe laser. Detailed derivations for the 
simplified cavity model can be found in subsections |B 1| and |B 2| 
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b. Probe interactions 



Subsystems of the qubit cavities corresponding to Z-probe interactions have the following limit model for the 
atom-cavity interaction: 

Q 12 = (Z< W ,0,0), Q 21 = (z^ ,0,0), Q 24 = [z^ , 0, 0) , Q32 = (^ (Q3) ,0,0). (C2) 
The analogous X-probe qubit cavities replace the substitution Z^ % > — > X^ l \ Detailed derivation of the probe 



interaction model is provided in subsection B 1 



Raman transitions 



The bit-flip network Raman systems are: 

Q 11 = (i,V7^ 1) 4An(w), Qis= (i^^AnW 1 
Q22 = ( i, yfrag* \ Unw) , q 23 =( 1, ^(f >, i Anw 



2) 



31 = (i,V7^? 3) ^An^), Q 33 = (i,Vi^ 3 \l^r Q3) ) 



(C3) 



where we note that the indices are connected slightly differently for Q\ than Q2 and Q3, following the diagram in 
Fig. jHJ The analogous phase-flip Raman interactions replace Ohr — > ~^( a hr + o-gr) and a gr — > -^(o~hr — o-gr) in the 
above. Details of the Raman transition model is provided in subsection |B 2| As explained in section [CM} the particular 
choice of coupling operators in each of the above helps to prevent Stark shift dephasing of the protected codeword 
without additional, Stark shift-compensating atomic levels. 



d. Coherent displacements 

Several vacuum inputs to the network will be displaced with Weyl operators, representing constant amplitude, 
coherent state inputs. For a given displacement, d, such inputs are represented with: 



W d = (I,d,0) 



(C4) 



e. Beamsplitters 

All beamsplitters in the network are taken to be 50/50, mixing the inputs of two fields: 

B = 



1 

x/2\ 



1 



\/2 \f2 



(C5) 



/. Relays 

Relays are modeled by a simplified four-state, four mode model. In the network there are two relays R\ and i? 2 , 
each of which has the decomposition Rk = Rki EB Rk2- The state of the relay is controlled by the field inputs to Rk2, 
while the state of the relay controls the routing of the field inputs to Rui, hence 




Tj(Rk) 
h 


rr(Rk) 


,0,0^ 


- „{Rk) 

(Rk) 


(Rk) ' 

Yf(Rk) 
h 


,0,oj 



(C6) 

for k = 1, 2, where the superscripts keep track of operators on different component states. Detailed physical modelling 
for such quantum relays and derivation of the simplified scattering model can be found in [5]. 
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g. Bit- /phase- flip errors 
The single-qubit error processes that the network seeks to correct may be described with 

= (l,VTx^\o) 
for a bit-flip error process with mean rate F on qubit i, and 

4 Qi) = (i,VrzW\o) 

for a phase-flip error process on qubit i, also with mean rate F. 



(C7) 



(C8) 



Network calculations 



The main task is to calculate the 'correction' network G ffl G' , with the subnets G and G' defined as follows. 
Parts of G : 

G p = R 12 < B < ((Q 12 < Q 21 ) ffl (7, 0, 0)) < B < (W^ a ffl (7, 0, 0)) , 

G f = (On ffl Q 3 i EH Q 22 ) < (S ffl 2 (7,0,0)) < (i? n ffl (7,0,0)) < (W^ ffl (7 2 , 0,0)) , (C9) 
representing the 'probe' and 'feedback' network subsystems of G — G p ffl Gf, respectively. Correspondences in G' : 

R22 «• R12, R21 «• Tin, Q23 «• Qu, Q33 «• Q22, Q13 Q31, Q24 «■ Q21, Q32 Ql2- (C10) 
The overall structure of the network is (with a trivial rearrangement of field indices) 



N 



diag(S( p \S'( p \s( f ),S'W), 



V 



7» 
L'(p) 

L(/) 

£/(/) 



\ 



H (p) + H'(p) + H« + H'« 



(Cll) 



/ 



In this calculation, the network components representing the bit- or phase-flips errors, E^ 1 ^, may be simply concate- 
nated to the 'correction' network N immediately before the derivation of the full system's master equation. 



a. Probe network 



Plugging in component models into G p and G' we arrive at 



G n 



G' 



(12) n (fll) + £ (12 )(7 gl) £ (12) n (Hl) + 0(12)^1) 

_£(i2) n (fli) _ Ods)^ 1 ) _ (i2) n f 1} - E^a { g f ] 

(32 )n (fi2) + E (32 )(J (R2) £(32 )n (*2) + (32) a g2) 

_ E (S2) U (m _ 0(32)^2) _ (32) n (fi2) _ £(32)^*2) 



V2 



n { g R1) 0^+a[f ] E^ 
-a g R1) 0^ -U[ R1) E^ 

nf 2) o( 32 ) + <7g 2) 7;( 32 ) 
-ag 2) o( 32 ) - nf 2) £( 32 ) 



(C12) 



where, for the bit-flip correcting network, E^ = Z^Z^ + 1 and OW) = Z<WZ<*'' - 1 arc proportional to 
two-qubit parity projectors. For the phase flip network, E^ = + 1 and = X^X^) _ \ 
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b. Feedback network 



Plugging in the bit-flip component models into G/ and G^ we arrive at 
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/— - (Q2 
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A sp 3 rrW*) 

2 Z^i=l 11 r 



/r (Q2)* n (Rl) 



A sp 3 nW*) 
2 Z^i=l 11 r 



Im 



V2< 



Im 



(Ql)* n (Rl)_ \ 



r (Q3)* n (Rl) 



9'' 



(Ql)* TT (R2) 



(Tor 



(Q2)* n (fl2) 



(7 



(Q3)* n (R2) 



\ 
/ 

(C13) 



We now perform the adiabatic elimination procedure on G/ffl Gy to isolate the Raman dynamics on the cavity 
qubits. To do this, we will exploit the trick described in section |B 3|in order to scale both the input coherent amplitude 



and excited state detuning, redefining j3 — > kj3, A — > fc 2 A, and G/ ffl G' f — > G^ ffl GV W ' with scaling parameter k 



Following the procedure in section B 3 wc construct the systems 
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/ 



Gf = 
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(fc)/ 



(Ql) n (fll) 

U gr ILg 



/9^(Q 3 ) n (Ri) _ rr (Q2) TT (fli) 

7Q1) , jx02) 
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P ) +2Im(/ffc/3(4? 1) *nf 1) -\ 



2Im 



(Q2)* n (R2) (Q3)* U (R2) 



) 



Define G 



limiting system 



G^ ffl &p' . Straightforward application of the usual limit theorem [13] taking k 



n( ffl )(nf) + n[ Q2) ) - 2ni m) n( Q3) - nf 2 )(nf) + 4 Q3) 



9n (R2) TT (Ql) 



(C14) 
oo yields the 

(C15) 



where the Raman interaction strength is now set by ft — 7|/3| 2 /2A. Again following the procedure of section 
construct Gp, the asymptotic approximation of G^p ffl Gf for sufficiently large values of k. 



B3 



we 



Gf = dia; 



^g R1) 
-^m ( g R1) 
-k(3u[ R2) 



,n( V2X^n g R x )nf 2) + x^n^u^ - V2x^u (R1) n^ 



g g 



h xi g 



^(nf) + nj Q2) ) - 2nf ^nf) - nf 2 )(nf 2 ) + n^ 3 )) - 2nf 2) ni Q1 



(C16) 



The analogous phase-flip network simply replaces with and 2II^ 4 ' ) — > /±jW'' in the Hamiltonian above. 
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TABLE I: AC Stark Shifts of error states AC Stark shifts for each three-qubit basis state in the bit-flip network described 
According to the Hamiltonian in equation (C18l, for joint each qubit state (Q1,Q2,Q3) and corresponding relay 



in section 



C3 



states that correctly reflect its parities (R1,R2), the AC Stark shift experienced by each qubit is listed under the columns (SSI, 
SS2, SS3). As designed, each joint qubit state will experience an aggregate — 2n AC Stark shift, as long as the error state is 
correctly represented in the state of the relays. 



Network master equation 



Concatenating the asymptotically approximate 'correction' subsystems calculated above with the error processes, 
we immediately arrive at a master equation for the bit-flip network (see section A 2 1 : 



dpt 
dt 



[H, Pt]+^\ L ]Pt L* - ^L*L jPt - ^PtL*Lj 1 



(C17) 



H 
Li 



. X (Q2) n (.Rl) n (fl2) 



n( ffl '(nf) + u[ Q2) ) - 2n{ R1) u^ - nf 2 )(nf + n 



(Q3)> 
h > 



a 



^(-ajfo^-nr^)), 



L 5 = VfX^\ L 6 = VfX^ 2 \ L 7 = VfX^ 3 \ 



The phase-flip network simply replaces XW*' O and 211^ 



(Qi 



(C18) 

/ ± A"(Q l ) in the above. Note that we have 



omitted the field coupling operators from in the above master equation. This is because these terms serve only 
to decohere the ground states of the relays, which have no effect on our network dynamics if the relays arc initiated 
in ground state eigenstates. The Hamiltonian dynamics principally serve to flip the state of the individual qubits, 
depending on the state of the relays (first three terms); L\ and Li drive the state of Rl to reflect the joint parity 
of Ql and Q2, and L3 and L4 drive R2 to reflect the parity of Q3 and Q2. The final four terms of the Hamiltonian 
are physically recognizable as AC Stark shifts that necessarily accompany the Raman interactions. Although not 
immediately obvious, the feedback was designed to minimize the detrimental effects of these shifts, as described in 
section IC41 



4. AC Stark shift compensation 



As discussed in sections B 2 and C 3 c our feedback network induces both Raman interactions that help to correct 
errors and AC Stark shifts that alter the energy of each qubit ground state. Unfortunately, these Stark shifts tend to 
break the energy degeneracy between joint qubit states, which consequently causes rapid dephasing in the encoded 
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FIG. 6: Performance with AC Stark shift compensation. Decay of fidelity, (vPolptl^o), for several values of the feedback 
parameter Q — /3 2 y/2A, a = fi/8, and T — 0.1 when extra, AC Stark shift compensating qubit levels are included. The fidelity 
decay of a single, bare qubit also suffering bitflip errors at rate r — .1 is also shown. 



codeword. There are two ways to counteract this. One is to design the network such that the net effect of these shifts 
is null and the other is to modify the components such that the AC Stark shift is nulled on a qubit-by-qubit basis. 

The network described in section |C 3| takes the former approach. As summarized in table |TJ by choosing the 
intensities of and the particular transitions targeted by the Raman beams, the feedback is designed such that the 
aggregate AC Stark shift experienced by any joint qubit state is equal to —20, as long as the relays accurately 
reflect the joint parities of the qubits. Thus, codeword dephasing occurs only in the interval between changes in 
the qubit parities and the accurate reflection of these changes in the relays. Consequently, the amount of Stark 
shift-induced dephasing may be decreased by increasing the probe amplitude a. This approach yields sub-optimal 
storage fidelities, as some codeword dephasing will occur for finite a, but has the advantage of significantly simpler 
components. Representative performance curves are shown in Fig. |4j 

The latter approach was discussed in some detail at the end of section |B 2| By presuming an extra pair of atomic 
excited states in each qubit with equal but opposite detuning from the Raman beams, the AC Stark shifts may be 
completely cancelled on a qubit-by-qubit basis. (Alternatively, an additional pair of feedback Raman beams with 
equal but opposite detunings would have an equivalent effect) The resulting network master equation would have an 
identical form as equations (C17) and (C18), but without the final four terms in the Hamiltonian. Such a network 



yields better storage fidelities than the network described in section [C 3| but requires significantly more complicated 
and precisely tuned component devices. Representative performance curves produced by this approach are depicted 
in Fig. |6j which were calculated for the same parameters as for Fig. [4] 
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